High harmonic spectra via dominant interaction Hamiltonians 
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We formulate the concept of dominant interaction Hamiltonians to obtain an integrable approxi- 
mation to the dynamics of an electron exposed to a strong laser held and an atomic potential leading 
to high harmonic generation. The concept relies on local information in phase space to switch be- 
tween the interactions. This information is provided by classical integrable trajectories from which 
we construct a semiclassical wave function. The high harmonic spectrum obtained is in excellent 
agreement with the accurate quantum spectrum. The separation in the atomic potential and laser 
coupling interactions should facilitate the calculation of high harmonic spectra in complex systems. 
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High Harmonic Generation (HHG) is one of the basic 
processes of non-linear light-matter interaction involving 
an electron under the simultaneous influence of a strong 
laser field and an atomic potential. Initial experimen- 
tal observations with atoms [l|, [|[ were soon followed by 
theoretical work 0, for an early review, see, e. g., 
More recently HHG has been also investigated in 
molecules [H, 0] and clusters [1, @ . 

The enormous impact of HHG up to recent proposals 
for imaging of molecular orbitals [lfj and the generation 
of attosecond pulses [ll[ is not the least due to a very 
simple description with the so-called three-step model 
[H, H| ■ According to this model the electron tunnels out 
of the combined nuclear plus laser potential, and is then 
accelerated and pushed back to the nucleus by the laser 
field where its energy can be converted upon recombina- 
tion into a high energy photon. This is possible up to 



a maximal cutoff energy of E n 



=3.17[/„ + /„, where 



Up = £ 2 /(4w 2 ) is the ponderomotive energy in terms of 
the maximal laser field strength £ and the photon fre- 
quency uj, while I p the ionization potential. E max was 
derived by considering classically the excursion of the 
electron in the laser field only [3|. On the other hand 
in purely classical HHG calculations the cutoff energy 
does not play a role since HHG is a quantum mechanical 
interference effect, essentially between electronic quan- 
tum amplitude in the ground state and in the contin- 
uum driven by the laser, as the three step model re- 
veals qualitatively. With a full semiclassical calculation 
based on trajectories in time, excellent agreement with 
the quantum result is achieved fl2l ]. However, this also 
implies that a quantitative description of HHG can only 
be achieved by solving numerically the non-separable dy- 
namics of an electron in the Coulomb field of a nucleus 
with potential V a and in the external laser field (Vs). 

In view of the universal relevance of HHG as men- 
tioned above it would be very desirable to have an ana- 
lytical approach which is also quantitative. To this end 
we introduce the idea of dominant interaction Hamilto- 
nians (DIH) to disentangle the non-separable dynamics 



by splitting it into spatial regions where either one of 
the two potentials dominates which is then taken as the 
only interaction in that region. For regions where the 
laser dominates (far away from the nucleus) the dynam- 
ics is simply that of a free electron in the laser field, 
the so called Volkov-dynamics, governed by the Hamilto- 
nian = p 2 /2 + Vg. In the opposite case we have the 
electron only under the influence of the Coulomb poten- 
tial V a which is trivially integrable for one electron with 
H a =p 2 /2 + V a . 

To demonstrate how the concept of DIH works, we will 
use the simplest realization of HHG, a one dimensional 
electron dynamics in a soft core potential V a along the 
linearly polarized laser field, defined by the Hamiltonian 
(atomic units (a.u.) are used unless stated otherwise) 



H = |- + V a + Ve, 



(1) 



where V a (x) = — (x 2 + a) -1 / 2 with a = 2 a.u. such 
that the ground state energy agrees with hydrogen 
(— Eb = 1/2 = I p ). The laser interaction is defined as 
Vs = x£ cos wi, where for convenience we use here a 3.5- 
cycle laser pulse with £ =0.1 a.u. and w = 0.0378 a.u. 
[12| . Our observable of interest, the HH spectrum 
a(u) can formulated as the Fourier transform c(lo) ~ 
J dt d(t) exp[iuji\ of the dipole acceleration 



d(t) = -(*(t)\dv a /dx\*(t)) 



(2) 



The exact ^f(t) is the solution of the time-dependent 
Schrodingcr equation with Hamiltonian Eq. ([1]) and leads 
to the familiar HH spectrum shown in Fig. [1] (a) which 
was extracted from the acceleration in Fig. [2] (a). It 
has been obtained under scattering conditions, i.e., with 
a Gaussian electron wavepacket ^i(x, 0) = (x\g(\Vi)}, 
where 



(x\g(wi)) = ( - 



exp 



--^ix-qif+^Piix-qi) 



(3) 
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FIG. 1: (color online) HH spectrum calculated from Fourier 
transformation of the dipole acceleration, (a) full quantum 
result, (b) semiclassical result for the full potential with 10 6 
trajectories, (c) DIH result with dominant interactions with 
10 4 trajectories. The cutoff for scattering from an ion under a 
laser field, -Ecutoft = 2U P +I P , determines the highest harmonic 
order, JV max = [-^cutoff/wj = 105. It is shown with a vertical 
(red) line. 



We have introduced the short notation w = (p, q) for 
a point in phase space. The wavepacket ^ with width 
parameter 7 = 0.05 a.u.. is initially located at rest (pi = 
0) at a distance of q. L = £/uj 2 (70 a.u.) from the proton 
which corresponds to half the quiver amplitude. 

Next, we construct the semiclassical HH spectrum. 
While the full quantum spectrum only serves as a refer- 
ence for accuracy of our approximation, the semiclassical 
propagation is part of the DIH approach to be developed 
since the latter requires local information in phase space 
as we will see. Such information is contained in the clas- 
sical trajectories underlying the semiclassical propagator 
developed by Herman and Kluk [l3| . (see also [3] and 
0), 



K(x,x',t) 



d 2 w 
2tt 



(x\g(w t ))VRe lb (ff(wo)la-') , (4) 



with Gaussians (x\g(w)) which have for convenience the 
same width as the initial state Eq. ([3]). The interpreta- 
tion of Eq. is straight forward: The quantum transi- 
tion amplitude from point x' at time t = to x at time 



FIG. 2: Dipole acceleration d(t), from which the HH spectrum 
shown in Fig. [1] has been obtained. 



t is constructed through classical trajectories which start 
at wo at time t = and reach under the dynamics of 
the Hamiltonian the phase space point w t at time t. The 
preexponential weight factor of such a trajectory in phase 
space is given by 



R 



-det 
2 



l PtPa 



l qtqo l l m qtPo > m vtqo 
7 



(5) 



which is composed out of the four blocks m a \y = da/db 
of the monodromy matrix [l3| . Note that for the present 
one-dimensional case, the m a b are scalars and no deter- 
minant has to be taken. The semiclassical amplitude is 
then given by "\/i?e lS ', where S(t) is the action along the 
trajectory. The integration is performed over all phase 
space points Wo which serve as initial conditions of clas- 
sical trajectories w t = (p t = p{p n , qo,t),q t = q(po, Qo, *))• 
Convergence is achieved with a finite number of trajec- 
tories through the Gaussian envelope, as illustrated in 
Fig. [3] The HH spectrum obtained with the semiclassical 
wavefunction ip(x,t) = J dxK (x, x', t)ipi(x') is in excel- 
lent agreement with the quantum spectrum, see Fig. [1] 
While providing a lot of insight into the dynamics which 
creates HH the semiclassical approach is at least as nu- 
merically involved as solving the Schrodingcr equation 
directly, since the trajectories w t cannot be obtained an- 
alytically, and moreover, the full classical dynamics of 
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FIG. 4: Soft Coulomb V a and laser potential — x£. Shaded 
areas represent the switching regions. 



FIG. 3: Distribution of initial conditions (dots) and analytical 
conditions for switching (lines, see Eq. ©). 



this problem is chaotic [12|, although the HH spectrum 
is very regular. 

This underscores the motivation for the DIH concept, 
where for each dominant interaction, the trajectories are 
ideally known analytically, or at least can be obtained 
with little numerical effort. The key point of the DIH ap- 
proach is to define an appropriate phase space boundary 
between the electron dynamics governed by the atomic 
potential V a and the laser potential Vg. For the present 
one dimensional case, the boundary reduces to isolated 
points w c in phase space. The physical process we have 
to describe is the trapping of the freely oscillating elec- 
tron in the laser field, due to the atomic potential. This 
will be most likely if the electron is slow close to the nu- 
cleus, i.e., we set p c = 0. The trapping region is defined 
in a natural way as the range where the soft-core poten- 
tial is always stronger than the laser potential. For V a 
and Vs from Eq. ([I]) , this region results to be the interval 



with x c = 3.0083 (see Fig. SJ). 



We are now in a position to generate the HH spectrum 
with the DIH approach. Although this can be done com- 
pletely analytically [lH, we prefer to use here the same 
propagation scheme of classical trajectories as used for 
the semiclassical HH spectrum presented above. This al- 
lows for strict comparison of the full semiclassical (and 
quantum) spectrum and the one to be calculated with 
DIH. The DIH result agrees remarkably well with the 
exact quantum spectrum as can be seen in Fig. [T] and 
Fig. El. 

In order to understand why the DIH approach works 
so well, it is instructive to analyze the initial conditions 
of the trajectories which switch and get trapped. A little 
thought reveals that they form bands in the initial phase 
space (see Fig. [3|). To see that we recall that the trajec- 
tory for an electron in a laser field with initial conditions 
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FIG. 5: Initial conditions for trajectories: switched trajecto- 
ries (diamonds), trapped or stranded trajectories in the full 
potential case (filled circles) and analytical conditions (lines 
in the inset, see Eq. ([9])). 



q(0) = qo and p(0) = po reads 

p(t) = po sin (uit) , 

q(t) = q + p Q t + — [cos (uit) - 1] 



(6) 
(7) 



From the condition p(t c ) = follows that t c sa nn/uj with 
n = 1, 2, . . . . Then, the conditions for a switch from Vg 
to V a are 

nir £ 

1c =Q0 H PO + COS (77,71-) - 1 

Ul Ul £ 



nir 



-qo 



-Po 



1 - (-1)" 2£ 



(8) 



This implies that the initial phase space points are given 
by lines p^do) with a width Apg = w/(nir) as illus- 
trated in Figs. [3] and [5] The explicit expression follows 
from rearranging Eq. ((8|), 

(n), x w ( 1 -(-!)" 2£\ 

Po M = qo 75 2 • 9 ) 

ri7r V 2 ui A 



4 



Interestingly, for the full interaction, the initial con- 
ditions for those trajectories which get trapped lie on 
the same phase space stripes (Fig. [5]). We may conclude 
that the DIH switching condition describes the dynam- 
ics relevant for HHG quite well. The small differences 
in the initial conditions can be attributed to the (small) 
attraction by the nucleus which the electrons from the 
full classical trajectories feel on the way inward. 

This observation also explains the long standing puz- 
zle why the Coulomb long range nature of the potential 
plays only a minor role: In fact, HH spectra in qualita- 
tive agreement with experiments have been also calcu- 
lated with zero-range potentials V a oc S(x) [17j . In the 
DIH approach, the switching condition would change to 
x c = but the structure of the initial manifold leading 
to switching remains the same. 

To summarize, we have introduced the concept of dom- 
inant interaction Hamiltonians (DIH) to simplify the the- 
oretical description of high harmonic generation by split- 
ting the problem into two integrable ones: the electron 
under the influence of the laser field and the electron un- 
der the influence of the atomic potential. We construct 
the HH spectrum semiclassically by using classical tra- 
jectories: They feel the force of the laser or of the atomic 
potential and the force is switched at the phase space 
boundary defining the dominance of each of the two in- 
teractions. The dynamics is integrable under either of 



the two interactions reducing greatly its complexity with- 
out loss of accuracy of the spectrum. The simplification 
manifests itself in the fact that up to one million tra- 
jectories are necessary in the present example of HHG to 
converge the spectrum fully semiclassically while a factor 
of 100 less is sufficient to converge the DIH spectrum. 

Moreover, the DIH provides a natural dynamical ex- 
tension of the simple man's approach: In the latter the 
wave functions of the electron in the laser field and the 
electron in the ground state of the atomic potential only 
are simply coherently added to produce qualitatively the 
HHG spectrum [l8[ . Here, with the help of DIH we have 
provided a framework how to dynamically populate one 
of these states (the ground state) while starting initially 
with the other one. 

While the HH spectrum presented here can be ob- 
tained fully analytically [l(| , the main thrust of the DIH 
lies in the perspective to describe HHG in more com- 
plex systems with the simplification of separating the in- 
teractions. Presently, one has to resort to semiclassical 
techniques since the switching condition is local in phase 
space. We will explore the possibility of a corresponding 
quantum condition in future work. 
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